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Abstract. Silk damping at redshifts 1.5 x 10'* < z < 2 x 10^ erases CMB anisotropics on scales 
corresponding to the comoving wavenumbers ^ < k < 10^ Mpc"^ This dissipated energy is gained 
by the CMB monopole, creating distortions from a blackbody in the CMB spectrum of the ju-type and 
the /-type. We study, using Fisher matrices, the constraints we can get from measurements of these 
spectral distortions on the primordial power spectrum from future experiments such as Pixie, and how 
these constraints change as we change the frequency resolution and the sensitivity of the experiment. 
We show that the additional information in the shape of the /-type distortions, in combination with 
the yU-type distortions, allows us to break the degeneracy between the amplitude and the spectral 
index of the power spectrum on these scales and leads to much tighter constraints. We quantify the 
information contained in both the yu-type distortions and the /-type distortions taking into account 
the partial degeneracy with the ^'-type distortions and the temperature of the blackbody part of the 
CMB. We also calculate the constraints possible on the primordial power spectrum when the spectral 
distortion information is combined with the CMB anisotropics measured by the WMAP, SPT, ACT 
and Planck experiments. 
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A Recipe for the calculation of CMB spectral distortions for a general energy injection 



1 The spectrum of CMB 

In the early Universe, at redshifts z » 2 x 10^, there is almost perfect thermal equilibrium between 
photons and electrons/baryons which maintains the spectrum of the cosmic microwave background 
(CMB) to be a blackbody spectrum even in the presence of enormous energy injection such as dur- 
ing electron-positron annihilation. This prediction of the standai^d big bang cosmological model 
was confirmed by the Far Infrared Absolute Spectrophotometer (FIRAS) instrument on the Cosmic 
Background Explorer satellite (COBE) [1] which found that the CMB is indeed blackbody to high 
precision. If there is any energy injection into (or cooling of) the primordial plasma at z > 2 x 10^, 
Compton scattering is able to very quickly redistribute this excess (or deficit) of energy over the entire 
spectrum of photons restoring the equilibrium Bose-Einstein distribution [2] with chemical potential 
parameter // and occupation number n{x) = l/(r'^^^ - 1), where x = hvlk-QT, h is Planck's con- 
stant, is Boltzmann's constant, y is the photon frequency and T is the temperature of photons and 
baryons. The chemical potential is in turn driven to zero by photon production [2] by bremsstrahlung 
and, more importantly for a low baryon density Universe such as ours, by double Compton scattering 
[3]. Electrons are also maintained at the equilibrium Maxwellian distribution by Compton scattering 
with the photons [4, 5] whose number density exceeds that of the electrons by a factor of ~ 10^. 
Coulomb collisions efficiently maintain equilibrium between electrons and ions in the entire redshift 
range of interest to us and they can be assumed to have the same temperature, defined by the photon 
spectrum[4, 5], Eq. A.l. 

At redshifts z < 2 x 10^, bremsstrahlung and double Compton scattering become inefficient 
in creating photons, however Compton scattering is still able to maintain kinetic equilibrium (Bose- 
Einstein spectrum) until z » 2 x 10^. This, therefore, defines the era where it is possible to create 
//-type distortions with the n par^ameter related to the fractional energy Q (= b£lEy, where tsE is the 
energy density going into the spectral distortions and Ey = a^T^ is the energy density of radiation 
and fiiR is the radiation constant) injected into the radiation by [2, 6] /i = lAQ. To calculate the 
jU-type distortions it is necessary to calculate precisely the suppression of ju behind the blackbody 
surface at z « 2 x 10^. This is of course possible with the numerical codes such as CosmoTherm^ 
[7] and KYPRIX[8], the former code includes precise calculation of distortions arising from the 
energy injection due to Silk damping. Sunyaev and Zeldovich [2] found an analytic solution for 
the suppression factor or blackbody visibility e"'"^^^ and this solution was recently made even more 
precise [9], allowing a fast and accurate computation of //-type distortions. These analytic solutions 
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and the recipe for using them are given in Appendix A. COBE 95% confidence level limit on the jj. 
parameter is < 9 x 10^^ [1]. 

At redshifts z < 2 x 10^, Compton scattering is not sufficient to maintain a Bose-Einstein 
spectrum in the presence of energy injection and is able to move the spectrum only partially towards 
the equilibrium creating intermediate or /-type distortions at redshifts z> 1.5 x 10^. The distortions 
in this epoch must be calculated numerically by solving the Kompaneets equation [10]. This was 
done in Ref. [11] and the results are publicly available. The recipe for using these results to calculate 
the /-type part of spectral distortions for a general energy injection scenario is given in Appendix A. 
We refer to [11] for a detailed discussion of the /-type distortions and how they can help distinguish 
between different energy injection scenarios, for example energy injection rate which is exponential 
in redshift such as particle decay, and that which is a power law such as Silk damping. We should 
emphasize that an /-type component is inevitable for power law energy injection such as Silk damping 
(which itself is unavoidable) while a particle decay can happen entirely in the ji epoch leading to 
negligible /-type distortions. Thus the presence of absence of /-type distortions together with the 
shape of the /-type distortions is a powerful discriminant for different energy injection mechanisms. 

At redshifts z < 1.5 x 10^, comptonization is minimal, and the solution for the distortions is 
given by a 3'-type distortion [12], which are also created at lower redshifts when the CMB photons 
are scattered in the clusters of galaxies by hot electrons. The 3'-type distortions are expected to be 
dominated by low redshift contributions coming from the epoch of reionization, with y parameter 
(see Eqs. 3.1,A.8) given hy y ~ 10~', and warm hot intergalactic medium [11, 13-20] with y ~ 
10"^ - 10"^. The .y-type distortions are therefore hard to predict and must be fitted as a free parameter 
during the Fisher matrix analysis. COBE 95% limit on the y parameter is j < 1 .5 x 10"^ [1]. However, 
a future experiment such as Cosmic Origins Explorer (COrE) [21] may be able to detect compact 
sources of y-type distortion and thus help estimate and clean the low redshift average j-type distortion 
contribution. 

There are spectral distortions, other than y, yU and /-type, which are created in the CMB, for 
example, recombination lines from the epoch of recombination [22-26] and in CMB anisotropics 
from resonant scattering of the CMB in metals lines during reionization and later [27], which we 
will not discuss here. We refer to a recent review [28] for a more complete discussion of the various 
spectral distortions in the CMB. 

2 Spectral distortions from Silk damping: Observing 17 e-folds of inflation 

Photons diffusing through the primordial plasma erase perturbations on small scales [31-33] and this 
effect, known as Silk damping, has already been observed in the CMB anisotropics by the Atacama 
Cosmology Telescope[34] (ACT), the South Pole Telescope[35] (SPT), and now the Planck experi- 
ment [36] on scales up to « 0.2Mpc"^ which is also the damping scale at z w 1200, where the 
anisotropics are suppressed but not completely erased. The wavenumbers coiTcsponding to the pho- 
ton diffusion length are ku « lO'^Mpc"^ (or multipole { « 10*^) at z ^ 2 x 10'' and ku w BMpc"^ 
{{ K 10^) at z = 1.5 X 10^, and on these very small scales the anisotropics are almost completely 
erased from the CMB and are therefore unobservable in the CMB anisotropy power spectrum. The 
energy stored in the perturbations (or the sound waves in the primordial radiation pressure dominated 
plasma) on the dissipating scales, however, does not disappear but goes into the monopole spectrum 
creating y, ji and /-type distortions, see Fig. 1. This effect was estimated initially by Sunyaev and 
Zeldovich [2] and later by Daly [37] and Hu, Scott and Silk [38]. Recently, the energy dissipated 
in Silk damping and going into the spectral distortions was calculated precisely in [39], correcting 
previous calculations and also giving a cleai^ physical interpretation of the effect in terms of mixing of 
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Figure 1. Power which disappears from the anisotropics appears in the monopole as spectral distortions. CMB 
power spectrum calculated using CMBFAST [29] with the undamped spectrum estimated by enforcing tight 
coupling in the code during recombination and cutting off the Une of sight integral at z = 800. See also [30]. 



blackbodies [39, 40] ^. The calculations in [39] showed that photon diffusion just mixes blackbodies 
and the resulting distortion is a j-type distortion which can comptonize into /-type or ju-type distor- 
tion, depending on the redshift. We can write down the (fractional) dissipated energy (Q = AE/Ey) 
going into the spectral distortions as [39, 40] 

do d r k^dk V Y^' 9 d r ^^d^ ^ r„9 „9i 

where @i{k) are the spherical harmonic multipole moments of temperature anisotropics of the CIVIB, 



t is proper time and P]{k) = qj^^^P^; « 1.45P^, - (A^27r2/)t3)(/t/)fco)"^"^^5d"s/dinA:(in/:/io)^ 
the amplitude of comoving curvature perturbation is equivalent to in Wilkinson IVIicrowave 
Anisotropy Probe (WIVIAP) papers [43], and Ry = Pvlipv + Py) * 0-4, Pv is the initial neutrino energy 
density and py is the initial photon energy density, but after electron positron annihilation [44], k^ is 
the pivot point, is the spectral index and d?is/dln^ is its running [45, 46]. The last approximate 
equality is valid in the tight coupling approximation when the energy in the multipole moments t >2 
can be neglected. 

It is possible to evaluate the time derivative in Eq. 2.1 explicitly using the first order Boltzmann 
equation for photons, yielding an exact expression in terms of photon quadrupole, dipole and higher 
order multipoles and baryon peculiai" velocity which is valid at all times [39, 40]. But in the redshift 



^See [41] for a slightly different way of calculating //-type distortions and also [42]. 
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range of interest to us, z > 1.5 x 10"^, the tight coupling solutions are quite accurate [47, 48], yielding 
the following expression for the energy injection rate [39, 49], 



f «p;,«,V-«»', (2.2) 



dz 4 dz J (2;r)3 ^ 
where is the damping wavenumber given by [31-33, 50], 

?2 



_L. fdz^iili^f^.a, (2.3) 



where R = ph is the baryon energy density, crj is Thomson cross section, me is the mass of 
electron, is the number density of electrons, c is the speed of light, and H is the Hubble parameter. 
In general the Eqs. 2.2 and 2.3 must be evaluated numerically. However in the radiation dominated 
epoch, /? «; 1 and for a power law initial power spectmm, d?is/d Ink = 0, we can evaluate the integrals 
analytically yielding, 

dQ 3.25Af d{l/ku^) „^,,(n + 1 



dz ytp'"' dz 



Once we know the energy injection rate, it is easy to calculate the yu-type and /-type distortions 
using the method given in Appendix A. We should also mention that the spectral distortions from 
SiUc damping may also permit us to constrain the primordial local type non-gaussianity on these very 
small scales [51, 52]. 



3 Fisher forecast results for Pixie-like experiments 

We will now calculate the constraints we can put on the initial power spectrum using CMB spectral 
distortions. Although the spectral distortions are sensitive to cosmological parameters other than 
the power spectrum, as can be seen from the equations, these sensitivities are relatively milder and 
we know the other cosmological parameters to high accuracy from CMB anisotropy and other data 
[34, 35, 43, 53-62]. We will therefore fix all parameters, except for the power spectrum, to the 
following values for flat ACDM cosmology[43]: CMB temperature Tcmb = 2.725 K, baryon density 
parameter Ob - 0.0458, cold dark matter Qcdm - 0.229, Hubble constant ho = 0.702, helium fraction 
0.24, eff"ective number of neutrinos[63] 3.046. 

We can write down the specti"al distortion of CMB, A/y, which will be measured by Pixie [64] 

as 

My = tll+yll + /y™P'"^(«s,Af,d?is/dlnA:). (3.1) 

The first term is the uncertainty in the temperature of the blackbody part of the spectrum which is 
not known a priori and must be fit simultaneously with the spectral distortions, second term is the 
y-type distortion which has contributions from low redshift intergalactic medium and is therefore a 
free parameter. The last term is the /-type -i- //-type distortions from the dissipation of acoustic modes 
and is a function of the power spectrum which we parmeterize by the amplitude, the spectral index 
and its running. The fomiulae for different terms in Eq. 3.1 aie given in Appendix A. If there is new 
physics injecting energy, there will be additional terms added to the above distortion. For example, 
if there is decay of particles before z = 2 x 10^, it will create additional //-type distortions and a 
term ///(,' should be fitted with // as a free parameter. Adding additional parameters will of course 
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degrade the constraints on the primordial power spectrum and we will discuss it briefly below. We 
should also include cooling of CMB due to energy transfer to baryons which cool faster than radiation 
[7, 22, 23, 49], however, it depends only on cosmological parameters which we have assumed to be 
constant and therefore does not affect the Fisher matrix. The cooling must be included in a precise 
calculation using, for example, Markov chain Monte Carlo to explore the likelihood. 

If Ay is the spectral resolution of the experiment and 61{v) is the noise in each channel, the 
Fisher matrix is given by (e.g. [65-67]), 



where vj is the center frequency of each channel, 9afi e {t,y,ns,A^,dnfi/dlnk), and the sum is over 
all frequency channels which we take to be from 30 GHz to 500 GHz. The upper limit of the usable 
frequency range will depend on how well we can remove the foregrounds which become more prob- 
lematic at higher frequencies. In the analysis below we will assume that the foregrounds have been 
removed at required precision (~ 61) using the frequency channels greater than 500 GHz which is the 
plan for the Pixie experiment [64]. In the Pixie proposal Av = 15 GHz, with a total of 400 frequency 
channels extending to 6 THz and the hope is that the large number of channels at high frequencies 
would help nail down the foregrounds to the desired accuracy. 

It is extremely important to study the impact the foregrounds have on the ability of Pixie-like 
experiments to detect spectral distortion but which is beyond the scope of present work. This paper is 
only a first step in quantifying the information content and the detectability of the spectral distortions 
with respect to constraining the primordial power spectrum on scales 8 Mpc"^ < ^ < lO'* Mpc"^ 

The covariance matrix is just the inverse of the Fisher matrix, cova/j = [^~^] ^ If we are in- 
terested in only first n parameters of the parameter vector 9, we can marginalize over the rest of the 
parameters by writing the Fisher matrix as 



where the sub-matrix A spans over the parameters we are interested in. The marginalized Fisher 
matrix is then given by F = A-BC^^B^ . It can also be obtained (if the parameters are non-degenerate) 
by taking the rows and columns of the parameters we are interested in from the covariance matrix 
and inverting the resulting sub-matrix. If we want to fix some parameter at a particular value instead 
of marginalizing over it, we just eliminate the row and column for that parameter from the Fisher 
matrix. The marginalized v - cr ellipsoids are then given by A9^FA9 = x^{v) where A9 = 9 - e/fiduciai 
and the Fisher matrix is also evaluated at the fiducial values of the parameters. For two parameters, 



Let us first consider the constraints we can obtain from spectral distortions alone. For this 
we take the pivot point = 42 Mpc"\ which is approximately in the middle of the /-type dis- 
tortions epoch, and constrain the amplitude and the spectral index on small scales (without run- 
ning and large scale power spectrum constraints). We thus have the parameter vector 9 = (ns(^o = 
42 Mpc"'),Af()to ^ 42 Mpc"'),?,y). We take the fiducial model n^{ko = 42 Mpc"^) ^ 0.96,A^()to = 
42Mpc"i) - 1.61x10"'^ and marginalize over t and y. The resulting 68% contours are shown in Fig. 
2 for Pixie [64] spectral resolution 6v = 15GHz and sensitivity 61{v) = 5 x lO^^^Wm^^Sr^^Hz"^ 
We also show the contours obtained by increasing the spectral resolution to 1 GHz or/and sensitiv- 
ity to 61{v) = 10~^^Wm~^Sr~^Hz~^. We should emphasize that these are completely independent 
constraints on the power spectrum on scales 8 Mpc"' < ^ < 10"^ Mpc~\ where there are cuixently 




(3.2) 




(3.3) 



(2) = 6.17. 
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1-a, Pixie:(Av=15 [GHz], AI=5 [lO'^* W m"^Hz"'Sr"']) Av=lGHz, AI^IO'^*' W m'^Hz"'Sr"' 




nj,(k„=42 Mpc"') "5(^0=42 Mpc'') 

Figure 2. Marginalized Fisher matrix constraints on amplitude and spectral index of primordial power spec- 
trum defined with respect to the pivot point = 42 Mpc Left panel shows constraints for different spectral 
resolutions and sensitivities of Pixie-like experiments. The labels for different contours are resolution in units 
of GHz and 51 in units of Wm^^Sr 'Hz (Av, 51). Right panel demonstrates degeneracies between different 
parameters for resolution 1 GHz and sensitivity Wm "Sr 'Hz ' . The curve labeled y free is the normal contour 
we expect marginalized over y and t. If we ignore the information in the shape of the /-type distortions and 
consider only yu-type distortions, as in studies so far, we get the curve labeled 'ju only'. If we add an additional 
free parameter, ju, which may arise from new physics such as decay of particles, we get the 'ju-free' curve. The 
curve labeled 'y-fixed' is the one we get if we assume we can predict and fix the y-type distortions from low 
redshifts. 

absolutely no constraints. Thus, although they may look much weaker compared to the constraints 
from CMB anisotropics [43], they still represent significant improvement in our knowledge of the 
cosmological initial conditions. 

The right panel in Fig. 2 demonstrates the degeneracies between different types of distortions. 
The curve marked 'y free' is the same as the (1,1) curve in the left panel, marginalized over y and t. 
The curve labeled '// only, no /-type' is the one obtained if we ignore that the /-type distortions have a 
characteristic shape with infonnation about the spectral index and divide all the energy released into 
j-type (z < 5 X lO'*) and ju-type (z > 5 x 10^) distortions, as in previous studies (c.f. Fig. 17 in [39]). 
Clearly it is not possible to constrain two parameters with one observable, /i, in this case and our 
two parameters are completely degenerate. Inclusion of the /-type distortions gives the 'y free' curve 
and converts the straight line contours into an ellipse demonstrating the infonnation contained in the 
shape of the /-type distortions. A different way of seeing this is by making the ^ parameter free, for 
example, as a result of new physics such as decaying particles at z > 2 x 10^. Note that we already 
have y and / free, so the constraints in the '// free' contours are coming solely from the shape of the 
/-type distortions. The fact that we can get any constraints at all in this last case demonstrates that the 
/-type distortions ai^e not completely degenerate with and cannot be mimicked by a combination of y, 
ju distortions and / parameter. Finally the curve labeled 'y fixed' is obtained by eliminating the row 
and column corresponding to y parameter from the Fisher matrix, thus fixing y. Comparing with 'y 
free' curve, this tells us how the presence of stars and galaxies (responsible for reionization, WHIM 
which give low redshift j^-distortions) in the Universe limits our ability to measure the primordial 
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Fiducial model nj,=1.018, dnd/dlnk=-0.022 Fiducial model 11^=0.965, dnd/dlnk=0.0 

— I 1 1 1 1 1 I 1 1 1 1 r 




0.96 0.98 1 1.02 1.04 1.06 1.08 0.9 0.92 0.94 0.96 0.98 1 1.02 

n,(k,)=0.002 Mpc"') nj,(k,)=0.002 Mpc"') 

Figure 3. Constraints (68% confidence level) by combining WMAP+SPT+ACT [34, 35, 43, 53, 54] with 
spectral distortions from Pixie-like experiments for two different fiducial models. Labels are same as in Fig. 2. 
A fiducial model with zero running has more power on small scales and therefore gives tighter constraints. 



power spectrum with the CMB spectral distortions. 

The model considered above is perhaps the most general model we can hope to constrain with 
the spectral distortions. A very restrictive model on the other hand is a model with running spec- 
tral index, which applies to both the large scale anisotropics and the spectral distortion. For this 
model, we take the pivot point at k - 0.002 Mpc"^^ as in the WMAP papers. So we can use the 
Markov chains provided by the WMAP team and combine the CMB anisotropy data with the spec- 
tral distortions to see how the spectral distortions improve the constraints on the spectral index and 
its running. We use the Markov chain with the running spectral index and including ACT [34, 54] 
and SPT [35, 53] data provided by the WMAP team^ [43] and use CosmoMC [68] to extract the 
covariance matrix for ?is, d?is/dlnA;, from it. We use two fiducial models, one with best fit WMAP 
values for the running model Ms = 1.018,d?is/dlnA: = -0.022, lO'^A^ = 2.345 and a second one with 
?is - 0.965, d?is/d In ^ = 0, 10^ = 2.43 to investigate how the change in fiducial model affects 
the constraints. We marginalize over y, t and A^ and give the 1 - tr contours for n^,An^lA\nk for 
different spectral resolution and sensitivities for the Pixie-like experiments in Fig. 3. It is clear from 
this figure that spectral distortion detection would be able to improve the constraints in this simple 
but restrictive model. The constraints are sensitive to the fiducial model. This is expected since a 
fiducial model with negative running will give much smaller distortions compared to a model with 
zero running and the effect is amplified because of the long separation of scales between the pivot 
point ak - 0.002 Mpc"' and the relevant damping scales at ^ > 8 Mpc"'. 

Planck CMB experiment's cosmology results^ are now available [36]. To estimate how an 
improvement in the large scale constraints on the piimordial power spectrum affect the information 
coming from spectral distoitions, we use the covariance matrix for baseline ACDM -i- running model 
provided by the Planck team[36] which in addition uses polarization data from WMAP [69] and high 

^http : //lambda . gsf c. nasa . gov/product/map/ current/ 

'^Based on observations obtained with Planck (http : / /www . esa . int/Planck), an ESA science mission with instru- 
ments and contributions directly funded by ESA Member States, NASA, and Canada. 
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Figure 4. Constraints (68% confidence level) obtained by combining the Planck+SPT+ACT+WMAP polar- 
ization [36, 43, 70, 71] with spectral distortions from Pixie-like experiments. Labels are same as in Fig. 2. 
Planck mean fit parameters are used for the fiducial model in the left panel of the plot. The right panel uses a 
fiducial model with zero running leading to stronger spectral distortions and tighter constraints. 



I data from SPT [70] and ACT [71]. For combining Planck results with spectral distortions, we use 
the Planck mean fit parameters, Q.\)h^ - 0.02225, Qcdm^J^ = 0.1205, /j - 0.672 and fiducial model 
with«s = 0.955,d?is/dlnA: = -0.015, ln(10i%) = 3.12,/to = 0.05 Mpc"^. The results are plotted 
in Fig. 4. We again demonstrate the effect of fiducial model in the right panel of Fig 4 where the 
fiducial model is taken to be = 0.958, dn^/d In yt = 0,ln(10'°Af) = 3.1, Atq = 0.05 Mpc"^ but the 
same Planck covariance matrix and other parameters as the left panel. By comparing Figs. 3 and 4 it 
is clear that there is additional and complementary information coming from the spectral distoitions 
irrespective of the improvements in the large scale constraints from the CMB anisotropics. 



4 Concluding remarks 



We have calculated the Fisher matrix forecasts for the possible constraints on the primordial power 
spectrum using CMB spectral distortions arising from Silk damping. These results, of course, come 
with the usual caveats associated with the Fisher matrix analysis. In particular, the spectral dis- 
tortions are a non-linear function of the primordial power spectrum parameters and are most likely 
non-Gaussian. In addition, adiabatic cooling of baryons [7, 49] reduces the signal and determines 
the smallest positive distortions that can be detected. Our analysis does not take this into account, 
since the constant terms drop out in the derivatives used to calculate Fisher matrices. For the cosmo- 
logical parameters preferred by Planck, the adiabatic cooling dominates over the Silk damping for 
d^s/d In ^ < -0.08, for the total /-type -i- //-type distortion, which is ruled out at high significance and 
our Fisher matrix calculations are a good approximation within the allowed region. Our calculation 
represents the first step in quantifying the information stored in the spectral distortions, in particular 
the /-type distortions and paves the way for a more careful analysis using Markov chain Monte Carlo 
techniques in the future. Fisher matrix analysis, in particular, shows approximately how much im- 
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provement in constraints can be expected if the spectral resolution or the sensitivity of Pixie is made 
better and also shows how these constraints are sensitive to the choice of fiducial models. 

We have also shown that there is important information in the shape of the spectral distortions 
coming from the /-type distortions, which has so far been ignored in constraint calculations, although 
this information was available from the numerical computations of the spectral distortions. In pai^- 
ticular, the /-type distortions are very important in breaking the degeneracy between the amplitude 
and the spectral index of the primordial power spectrum on the scales S < k < 10^ Mpc^^ More 
important than the error bars on different power spectrum parameters is the fact that with spectral 
distortions we will be able to extend our knowledge of initial conditions to completely new scales 
separated by many orders of magnitude from the information available from CMB anisotropy and 
large scale structure data. In inflationary terms, this amounts to extending our view of inflation from 
~ 6 e-folds at present to ~ 17 e-folds, which might well be a significant fraction of the full inflationary 
epoch. 



A Recipe for the calculation of CMB spectral distortions for a general energy injec- 
tion scenario 



To calculate the spectral distortions arising from a given energy injection rate, dQ/dz, where Q = 
AE/Ey is the fractional energy injected into the CMB, we must solve Kompaneets equation[10] in- 
cluding photon production and absorption due to bremsstrahlung and double Compton scattering [2]. 
The problem is however complicated by the fact that the electron temperature enters the partial differ- 
ential equation describing the evolution of photon intensity ly, or equivalently the photon occupation 
number n{x) = l{2hv^)Iy, itself depends on the photon spectrum [4, 5] 



T ~ 



4 Jnx^dx 



(A.l) 



At high redshifts, when Compton scattering is able to establish Bose-Einstein spectrum, an 
analytic solution for the evolution of ju parameter accurate to ~ 1% can be found and is given by 
[2,9] 



\ dz dz 



-Til) 



(A.2) 



where 
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(A.4) 



5.916Cabr 



fldC = neoO-jC—— — gdc(.Xe)IdC 



afs^BO /^B^CMB^ ^^^/ 



C = 0.7768, /dc ~ 25.976, gdc = 1-005 and gbr = 2.99, all quantities with subscript zero ai^e 
evaluated at redshift z = 0, o-fs is the fine structure constant and Or is the total radiation energy 
density parameter with relativistic neutrinos. The above equations assume that there is injection of 
only energy. If there is also significant injection of photons (other than bremsstrahlung and double 
Compton) then an additional term, -2.404^ must be added in the brackets in Eq. A.2, where A'^ is 
the fractional change in the number density of photons, z^ax should be taken to be sufficiently behind 
blackbody surface at ZdC * 1.96 x 10^, Zmax = 5 x 10^ should be sufficient for most energy injection 
scenarios, is the boundary of transition from //-type to /-type distortions and is discussed below. 
We have also accounted for the cooling of radiation because of energy transfer to baryons which cool 
faster than radiation with the expansion of the Universe and which has a simple expression before the 
start of the recombination [7, 22, 23, 49], 



cooling 



3 feB(»H +»He + »e) 

2 a^T\\+z) ' 



(A.5) 



where nu and ?iHe are the number densities of hydrogen and helium nuclei. 

For the /-type distortions we must solve the Kompaneets equation numerically. It turns out that 
the Kompaneets equation can be cast entirely in terms of dimensionless variables using Vy as the time 
variable. 



yy(Z,Zm}) 



»-'<.ini 



dz 



m^c H{\ + z)' 



(A.6) 



Solving Kompaneets equation with the initial spectrum corresponding to a 3'-type distortion, it was 
found[ll] that the spectrum starts deviating from j-type distortion at 1% level at = 0.01 and is with 
1% of a yu-type distortion at jy - 2. We thus define the boundaries of the /-type epoch by yyiO, z^) = 2 
and yy{0,Zy) = 0.01. At z > « 2 x 10^ we have yu-type distortions and at z < Zy w 1.5 x 10^ we 
have y-type distortions. 

The total distortion, excluding y-type, is now given by. 



2hv^ 



m IdQ dQ'°°''^^\ dz ^ 



(A.7) 



where all terms are evaluated at redshift z\^^ related to yy . by jy^ - yy{0, z\^^, the sum is over values of 



inj' 



jy finely sampled between 0.01 and 2, njiyy.) is the intermediate spectrum obtained by evolving j-type 
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distortion with energy 2ref fromjy - to . with Kompaneets equation, = CVy,-^j -3'y,_j)/2. The 
/-type spectra, sampled at intervals dvy^ = 0.001 for yy < I and Syy^ = 0.01 for 1 < Vy < 10 for 
Gref = 4 X 10~^ are available at http://www.mpa-garching.mpg.de/~khatri/idistort.html 
along with a Mathematica code which implements the above recipe. The above formula simply 
calculates the energy injected in each small redshift interval and adds the appropriate distortion to the 
total. With the redshift/j^ bins defined above the accuracy of the final spectrum is ~ 1%. The most 
time consuming part of the above calculation is the calculation of yyiO,Zini) as a function of Zmj, but 
it can be stored and reused if the normal cosmological parameters are not changing significantly. 
Finally, 3'-type , yu-type and f-type occupation numbers are^ [2, 6, 12] 



xe 



{e^ - 1)2 

(ev-l)2\2.19 

xe'^ 

{e^ - If 



+ 1 

e'' - 1 



- 1 



(A.8) 
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